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Abstract 

' We pursue the analysis made in [I] on the arithmetic area enclosed by m closed Brownian 

, paths. We pay a particular attention to the random variable 5 , niin2j ... i „ m (m) which is the 

. arithmetic area of the set of points, also called winding sectors, enclosed ri\ times by path 1, 

ri2 times by path 2, . . . , n m times by path m. Various results are obtained in the asymptotic 
| limit m — > oo. A key observation is that, since the paths are independent, one can use in 

. the m paths case the SLE information, valid in the 1-path case, on the 0-winding sectors 

arithmetic area. 



^ ! 1 Introduction 

O-; 

_C In [T], the asymptotic behavior of the average arithmetic area enclosed by the external frontier 

of m independent closed Brownian planar paths, of same length t and starting from and ending 
at the same point, has been obtained using a path integral approach [2 [3J. 

In the one path case, the random variable of interest happens to be the arithmetic area S n 
of the n-winding sectors enclosed by the path, from which the total arithmetic area S = ^2 n S n 
can be computed. A n-winding sector is by definition a set of points enclosed n times by the 
path. Path integral technics [3J give (S n ) = t/(2nn 2 ) but end up being somehow less adapted 
for (So), the average arithmetic area of the 0-winding sectors inside the path, i.e. of the set 
of points enclosed an equal number of times clockwise and anti-clockwise by the path. Indeed 
path integral cannot distinguish 0-winding sectors inside the path from the outside of the path, 
CN ' which is also 0-winding. Other technics have to be used, in the case at hand SLE technics [3], 

to get (S) = irt/5 from which (So) = 7rt/30 can be derived. It means that (So) = q(S — So) 
with q = 1/5. 

For m independant paths, the same path integral technics used in the one path case have 
now to focus [lj on the random variable S n (m), the arithmetic area of the n-winding sectors 
^ 1 enclosed by the m paths, from which the total area S(m) = J2n^n(m) can in principle be 

computed. A n-winding sector is again defined as a set of points enclosed n times by the m 
paths as illustrated in Figure [1] in the m = 2 paths case. One has found that the leading 
asymptotic term scales like lnm, namely that (S(m) — So(m)) ~ ^lnm with, as already 
stressed, no information on (So(m)) inside the m paths. 

In [5] on the other hand, the asymptotic behavior of the arithmetic area enclosed by the 
convex envelop of m closed paths was found to be (<S(n?,)) convex ~ y lnm, i.e. the same scaling 
as (S(m) — So(m)). One then concluded [T] that, necessarily, (S(m)) ~ y lnm and that 
(So(m)) is subleading. 

This particular scaling might be of interest for spatial ecology considerations where one 
asks about the increase in size of a natural reserve with the animal population (assuming 
of course that food is homogeneously and sufficiently available). By identifying an animal 
looking for food to a random walker one thereby obtains the logm scaling rather than the 
naive geometrical *Jm scaling. 
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Figure 1: On a square lattice two independent closed walks 1 and 2 with 38 steps each, starting 
from and returning at the origin. The winding sectors of each walk are labelled by their winding 
numbers {ni} (for walk 1) and {ri2} (for walk 2). The winding sectors enclosed by the external 
frontier of the two walks 1 + 2 are labelled i) by their joint winding numbers {n±,n2}, ii) by 
their total winding number {n\ + n2 = n}. The n = 0- winding sectors of interest, namely inside 
the external frontier of 1 + 2, correspond either to {0, 0} with at least one of the 0-winding 
sectors inside one of the paths, or to {ni, 11,2} with n\ + ri2 = and both n\ and n2 / 0. 
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In the present work we will give a more detailed analysis of the random variable S n {m). In 
particular we will pay attention to the random variable S fllB! ,„ nm (m), the arithmetic area of 
the set of points enclosed n\ times by path 1, rti times by path 2, . . . , n m times by path m, with 
Yli Tii = n. Happily enough, this variable can be tackled by path integral technics analogous 
to the one used in [IJ [3] provided that at least one of the n% ^ 0. Note, when J2i n % = 0> that 
Sn 1 ,n2,...,n m ( m ) contains a part of the inside 0-winding sector area So(m). To get a hand on the 
other part Soa—^O 71 )) with at least one of the rtj = 0-winding sectors inside the corresponding 
path i (otherwise one would be trivially outside the external frontier of the m paths), one key 
observation is that since the paths are independent one can use in the m paths case the SLE 
information (So) = q{S — So) valid in the one path case. 

From these considerations will follow (S(m)) and (So(m)). We will show in particular 
that, when m —> oo, the subleading (So(m)) remains finite. Also some information on the 
asymptotics of {S nijn2t o,...,o(m)) when both n\ and n<i — > oo will be obtained. Finally (So(2)) 
and the overlap between 2 paths, (2S(1) — S(2)), will be considered. This quantity might have 
some interest in polymer physics where polymers are modelised by Brownian paths. 

In view of unifying notations between [1] and [3], we always denote the n- winding sectors 
arithmetic area and total arithmetic area for m paths by S n (m) and S(m), which means that, 
from now on, S n (l) and S(l) stand for S n and S, the n-winding sectors arithmetic area and 
total arithmetic area in the one path case. 



2 Winding sectors 

2.1 Winding angle and propagator 

As stated in the Introduction, the arithmetic area of the winding sectors enclosed by planar 
Brownian paths can be obtained from their winding properties. Consider a path of length t, 
starting from and ending at and the angle 9 wound by the path around the origin O. The 
average of the random variable e tad over the set of such paths is 

(e^> = ^5 (!) 

where 

g a (f,r) = [ r 1 " r W(r)e"^o^W dT + io /o 9 M dT (2) 
Jr f (p)=f 

is the quantum propagator of a charged particle coupled to a vortex at location O. By sym- 
metry, it depends only on r 

1 2 +°° / r 2\ 
ft.Kf)= £ V«l T (3) 

k=—oo \ / 

where I\k- a \ 1S a modified Bessel function and 

Additional symmetry and periodicity considerations lead to Q a = Ga+i = Q\- a so that a can 

be restricted to < a < 1. We also set — = x so that 
~ ~ t 

+oo 

(e iad )=e- x J2 r lk _ a] (x) = G a (x) 

k=—oo 
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with Gq(x) = 1. Clearly areas are proportional to t, so we can, without loss of generality, set 
t = l. 

Let us rewrite G a {x) in a more suitable form : one has 

+00 

G a {x) = e~ x Y,h+ a (x) + {a -> 1 - a} (5) 

^ / +00 \ 
Observing that — ( e~ x ^ 4+a(«) ) = (J a _i(z) - we get 



fc=o 



= (J_ a (x) - 7 a (x) + I a -i(x) - h- a {x)) = e~ x (K a (x) + K X - tt {x)) 

(XX I 7T 

(6) 

Using the integral representation [6] of K a (x) 

1 r°° 

K a (x) = - du e~ xcoshu cosh(au) (7) 
J —00 

we deduce 

sin(avr) f°° 1 - e -*(i+c°sh u ) u 1 

G a (*) = — _ y_ dn i + coshu cosh 2 COSh((a " 2> } (8) 

Clearly 

lim = / du : — cosh — cosh((a )u) = 1 (9) 

x-+oc v 1 vr i_oo 1 + coshu 2 u 2' ' y ' 

so that 

„ , , sin(avr) Z" 00 e -^(i+cosh M ) 1 

1 _ G a (z) = — i — - / du cosh - cosh( (a - - u (10) 

7T j-00 1 + coshu 2 2 

Eq. (|10p will be extensively used in the following. 

2.2 Average area (S n (m)) of n- winding sectors 

Let us first consider the average arithmetic area (S n (m)) of the n-winding sectors labelled by 
their winding number n. For m brownian paths of same length unity, starting from and ending 
at the same point, one has [1] 

Z a (m) = ir H dx{l-{G a {x)) m ) = Y J {Sn{m))(l-e l2 ™ n ) (11) 
Rewriting G a {x) = 1 — (1 — G a (x)) leads, for n 7^ 0, to 



(S n (m)) 



£ da Z Q (m)cos(2vran) = ^j(-l)^ +1 (™^J £ da J™ dx (1 - G a (x)) j cos(2vran) 



where ( m ) is the binomial coefficient. Thus 



-■3 



(Sn(m)) = J^-iy +1 \Jjl n , j (13) 
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with 

"1 poo 



'0 JO 



I n j = — 7T J da J dx (1 — G a (x)) 3 cos(27ran) (14) 
Using {JTUJ) one rewrites I n j as 



dli,; \ $ r . 



- - (n 27r co ; h f j . + E / j cosh ui ^ 

where the ^Vj's follow from the integration over a 

• j odd = 2k + l 

" 2 2 ^ i 2 J J^J 1} ^ + k - n) {7r{2 N - I)) 2 + (ELi 

(16) 

• j even = 2k 

- gfrd>>«» (%^) ;g ( y »_.) _ + ; £Li _ 

(17) 

It is possible to compute I ni i exactly [J| 

^ f°° / 2JV-1 2iV + l \ 1 

n>1 ~ 2 y_oo l + coshu \(ir(2N- l)) 2 + u 2 (7r(2iV + l)) 2 + u 2 ) ~ 2vrn 2 [ ' 

On the other hand when n — > oo 

• j odd 

I n j ~ TIT C 7 ( 19 ) 



dx e- jx (K (x)) j (20) 

JO 

• j even 

(-i)(*+p j + 1)1 

/nJ ~ (2^ n> +2 dj (21) 

d,-=/ dxe-^iKoix))^ 1 du e~ xcoshu u t&nh- (22) 
Jo jo 2 

We finally obtain 

<s n (i)> = i nA 1 



2vrn 2 

2 2 3 1 

(Sn(2)) = W " /n ' 2 = ^°° W " 2^? d2 + ° ( ^ ; 
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and in the general case 



in f w m m \ t m \ t 

(S n (m)) = — 2 - I I n , 2 + I J n , 3 



n— >oo 



m 3 / m\ , / m\ \ , 1 . 



with the convention that (™) = if j > m, c?2 ~ 2.84 and C3 ~ 5.73. Clearly, and as discussed in 
the Introduction, we have no information so far on the 0-winding sector inside the m paths. To 
make some progress on this issue, we have to turn to (5 , nijn2i ... irim (m)), the average arithmetic 
area of the sectors enclosed n\ times by path 1, 112 times by path 2, . . . , and n m times by path 
m. 

2.3 Average area (S nitna) ... >nm (m)) 

Winding sectors can as well be labelled by the set {ni,ri2, • • • ,%} of the individual winding 
numbers rn enclosed by each path i. In line with section (|2.2p . one has 

poo I m \ 

Z auaa ,... tam (m)=Trj Q dx\l-Y[G ai (x)J = E '( S m,n 2 ,...,njm)) (l - e^tE^i 

(23) 

where in J2' the set n\ = n% = ... = n m = is excluded from the sum. So 



{£>m,n2,-,nm ( m )) 



jf dai... da m Z ai)Q . 2) ... )am (m) cos ^27r a i n ^ J 
/ dai... / da m Z ai>a2i ... )Q , m (m) TTcos(27raini) (24) 



where we have used Jq 1 da sin(7ra) sin(27rare) cosh na — = 0. Now {S ni ,n2,...,n m { m )) ls m ~ 
variant by permutation on the n^s, so one can focus without loss of generality on (S ni; ___ ;nj) o ; ___ ; o(m 
1 < j < m, with m, ...,rij 7^ 0. 

Rewriting G ai (x) = 1 — (1 — G ai (x)) in (|23l24j) leads to consider when m 7^ 

pi poo 

da i (l-G cti (x))co S (27ra i n i )= d Ui e^^+^^P^, m) (25) 
Jo J -co 

with 

uf + (1 - 4n 2 )vr 2 1 
(Ui ' " i} " K 2 + (vr(2n i + l))2)(u 2 + (vr(2n J -l))2) " 4^f (26) 

It follows that 

/•oo / p poo \ 

(S w ,..., % ,o,.,oH) = tt(-1) 3 ' ^ dx [] /^d^ e-^ 1+coshu ^ Pfam) (1 - Z(x)) m ~ J ' 

(27) 

where 

pi poo p— x(l+cosh«) 

/(x) = / da (1 - G a (x)) = / du (28) 

JO J —00 U + 7T 
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For example when m = 1, 2 one gets 



oo poo 



1 



<S ni (l)> = -tt/ dx I due-^ 1+cosh ^P(u,ni) 



u J-oo #kti\ 



OO /"OO 



{■Sm.,,(2). - <t/ / — ^ ld ^ 2 - P( Ul , ni )P(u 2 ,n 2 ) 
' 1 2 + cosh u\ + cosh u 2 



OO J — OO 



1 roo rco A*. A v 

2ixn\ J-oo J-oo 2 + cosh u\ + cosh ^2 

When nj — t- oo one obtains 

(S m ,... 1% -,o,..,o(m)) ~ 2i7r 2i-i n 2_ n 2 c i,™ 

where the constant 

/■OO 

dx e"^ (K (x))^ (1 - /(x)) m ^ 



o 



has to be evaluated numerically to the exception of 

C2 , 2 = |c(2) - 1 + |c(2) E na - - E na - ^i) 2 
z z fe=ii=i /z fc=ii=i + i 

3 Arithmetic area enclosed by m Brownian paths 

We are now in position to compute {S(m)), the average arithmetic area enclosed by the m 
paths. Obviously 

{S{m)) = £ '(S nu ...,n m (m)) + (5o,...,o(m)) (29) 

where (So,...^^)) is the area of the finite {m = 0, ri2 = 0, . . . , n m = 0} winding sectors inside, 
at least, one of the pathfl From (J23|), we get 

'(Sn 1 ,...,n m {m)) = f dai... / da m Z aii ... i£tm (m) (30) 

JO JO 

so that, using (f28|) . 

/•OO 

E '(^.....njm)) = vr / dx (1 - (1 - /(x)) m ) (31) 

JO 

In order to find (So,...,o( m )) first consider 

("OO 

E (Sn u o,...,o(m)) = vr / dx /(x) (1 - /(x))" 1 " 1 
■ ,~ Jo 



ni^O 



and, the sum 



A l = E (( S n 1 ,0,...,o( m )) + ( S 0,n 1 ,0..,o(m)) + - + (<So,...,0,rH ( m )>) 



m 



7T / dx/(x)(l-/(x)) 
1 / JO 



m—1 



'Remember that (<So(m)} = (So,...,o{m)) + X) (^ n i> ...,n m (jn)), with rij = and at least one of the m / 0. 
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where one has taken into account permutation invariance. Similarly one can consider 

A 2 = (( S n 1 ,n 2 A-A m )) + - + (So...,0,n 1 ,0...,0,n2,0...o( m )) + - + (£o,.,0,"i,n 2 ( m ))) 

ni,n 2 ^0 



V 2 

and in general 



dxf(xY(i-f(x)y 







A rn ^ (S ni ,n 2 ,...,n m (rn)) = r% dxf(x) m (l-f(x)) 



ni,ri2. ..,n m ^0 

with, obviously, 



771 POO 

£^ = tt / (1 - (1 - /(x)D = £ '(Sn u n 2 ,...,n m (m)) (32) 

i=l 

We are interested in (So,...,o( m )) : in the case of one closed path one knows from SLE [3] that 

„ <So(l)> 1 



En#0<5n(l)) 5 



(33) 



This means that for any point inside the path, with winding number n, q is the ratio of the 
probability to have a 0-winding to the probability to have an / 0-winding. 

Ai counts the points with only one non-zero winding number. It follows that the corre- 
sponding contribution to (<Sb,...,o(?Ti)) is necessarily q A\. Similarly, A2 counts the points with 
only two non-zero winding numbers. Since the m paths are independent, it follows that the 
corresponding contribution to (Sq ... o( m )) is Q 2 A 2- This line of reasoning generalizes to A^ : 
the contribution to (5 , o i ... i o(m<)) is q k A^. Finally 

771 

<S ,...,o(m)> = £ q% = $ (m) - %(m) (34) 
i=i 

with 00 / 

%(m) = 7r | o °° dx(l - (1 - (1 - g)/(x)) m J (35) 

Clearly YaLi A i m (j32[) coincides with <l>o(m). It follows that the average arithmetic area 
enclosed by the m paths is 

m 

(S(m)) = 2(1 + q l )Ai = 2$ (m) - * ff (m) (36) 
i=l 

When m — )• 00 (see eq. (|4"6j) in the appendix) one obtains 



(S(m)) = \ Inm - \ Inlnm + f (m ^ + Cj + o{-±=) (37) 

where C is the Euler constant. In Figure [21 numerical simulations for {S{m)) show that the 
agreement with eq. ([37|) is quite correct, even for not large values of m. 



•2 
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Figure 2: The average arithmetic area (5(m)): the crosses are numerical simulations (10000 
events) of closed random walks (10 6 steps for each one) on the 2D square lattice; m = 
4, 8, 16, 1024; the line is the analytical result eq. (|37p . The agreement is quite correct as 
soon as m > 16. 



Moreover, the asymptotic [1| of (S(m) — 5o(m)) when m — > oo is known to be 



Comparing Eqs. (i37l) and (|38l) . one deduces that the subleading 0- winding sector average arith- 
metic area 

lim (5 (m)) = - In- (39) 

remains finite in the m — > oo limit. 

Finally, as another illustration of the path integral formalism, consider the average overlap 
(25(1) — S(2)) of the arithmetic areas of two paths and (So (2)), the average 0-winding sectors 
arithmetic area of two paths. One has 



(S(m) — So(m)) = — lnm 




lnlnm + - (-\nV^r + c) +o(-==) (38) 




so that 



r 

vr(2-(l-g) 2 ) / 
Jo 



(25(1) - 5(2)) 



dx f(xf 




(40) 



Numerically 



(25(1) -5(2)) 
(S(l)) 



0.286 
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is close to what one would obtain if the paths attached in O were two circles of radius R: the 
overlap in unit of ttR 2 would then be 



1 2 

_ ~ 0.297 

2 7T Z 

Also, as far as (5*0(2)) is concerned, rewrite 

(S (2)> = (S(2)) - £<S„(2)> = 2(5(1)) - (25(1) - 5(2)) - ]T(5 n (2)) 



From Section E2] one has 



£(5 n (2)) = --7r / / dmdn 

n^O J J-ooJ-oo 



°° f°° , tanh ^ + tanh ^ 



' (2 + coshui + coshu 2 )(ui + «2)((2vr) 2 + (ui + u 2 ) 2 ) 



7f 

so that, using (5o(l)) = — and eq. (|4"0"|) . one finds (25o(l) — 5o(2)) to be 

30 



7T 



00 f°° duidii 2 / f§ tanh ^ + tanh ^ 



3 (2 + cosh m + coshu 2 ) Uvr 2 + n 2 )(vr 2 + u\) (m + n 2 )((2vr) 2 + ( Ul + u 2 ) 2 ) J 

(41) 

If the paths were not overlapping, then one would necessarily have (25o(l) — 5o(2)) = 0. 
The non vanishing result (14ip clearly indicates that the two paths do, on average, overlap as 
already seen in (|40]l . Note that the 0-winding sectors "overlap" (25o(l) — 5o(2)) is different 
in nature from the arithmetic area overlap (25(1) — 5(2)): the latter is a purely geometric 
overlap, whereas the former is more subtle since the superposition of the two paths destroy 
some original 0-winding sectors in each of the paths and create new 0-windings sectors for the 
two paths. 



4 Conclusion 

Path integral technics have been extensively used to tackle the issue of n-winding sectors 
arithmetic area of m Brownian paths. Some information stemming from SLE technics valid 
only in the one path case have also proved useful in the m paths case, merely because the 
paths are independent. 

Eqs. (|37p and (|39p are the main results of this paper. In particular the subleading 0- 
winding arithmetic area has been shown to remain finite in the asymtotic limit. Numerical 
simulations have nicely confirmed these asymptotics results. The overlap between two paths is 
also computed numerically. Applications to polymer physics will be studied in a forthcoming 
publication. 

5 Appendix 

In the appendix we derive the m — > oo asymptotic limit of & q (m) defined in (|35p 

poo 

<S> q (m)=TT dx(l-(l-(l-q)f(x)r) 
Jo 

with f{x) given in (|28p 

/oo g— x(l+coshu) 
-oo 7T + U Z 
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. ylnm 
betting x = — - — , 9q{m) becomes 



i » Trlnm f 00 , / . . . „i/lnra,,, 

*,N = — g~ ^ dy (^1 - (1 - (1 - g)/(iL_))r 

dy (l - e -» (42) 



7rlnm 



Using /(x) ~ — =-y — when x — )• oo, the integrand in (132" 

7T Z V X 



behaves, in the m — > oo limit, like 

if y > 1 

1 if y < 1 

i / \ vrlnm 

and so, at leading order, <P q (m) = — - — . 

1 

Focusing now on the subleading correction at order , namely 

ylnm 



^ m 1 '°°dy fl-e-^^V^ _ / V (1 -^V^ I (43) 



2 Ui \ J Jo 

let us first consider the y > 1 integration. One has to compute 



7rlnm f°° m 1 y \ Att 7rlnm f°° m 1 y j Air 

dy(l-g) — ~ — n / dy(l-g)- 



7r 2 y ylnm 2 J\ tt 2 V lnm 

where one has used that, because of the m l ~ y factor, y is peaked to 1 when m — > oo. One 
obtains 

a = ^= + o(-=) (44) 
V tt him Vmm 

Considering next the y < 1 integration one has to compute 



b 



7rlnm 


/ dye 
Jo 




2 , 




7rlnm 


/ dye" 
Jo 




2 , 




7rlnm 


/ dye" 


■a'm 1 - y 



4tt 
y In m 



/ 47T 

In m 



with a' — ' 



^ Jo 

dz 

Setting a'm}~ y = z with — dylnm = — leads to 

z 



ir 2 V lnm 



IT [ am _ z tiz TT I a'-m , /" a '"\ 
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At order , one has [In z e z ]^, m ~ — In a' (1 — a') and f£, m In z e z dz ps / °° In z e z dz — 

Jq Inz e~ z dz ps C — J a In 2 (1 — z)dz where C is the Euler constant. The last remaining 
integral is straightforward and finally 

bra I (W + C-a') (45) 

7T 1 — (7 7T 

Noticing that a'— = ^^^^ a, we are left with a + b rj — (In a' + C) so that 
2 virlnm 2 

$„(m) = -lnm - ^lnlnm + - ( ln(l - g) + In \/4t + C ] + o( 1 ) (46) 
2 4 2 y V 7T 3 y Vlnm 
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